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ABSTRACT 

Ocean  currents  can  cause  significant  and  interesting  effects  on  the 
intensity  of  underwater  sound  transmissions.  We  study  this  phenomenon  via 
the  parabolic  approximation,  beginning  with  conservation  laws,  and  derive  a 
family  of  equations,  each  of  which  is  valid  for  different  magnitudes  of  cur¬ 
rent  speed,  current  gradient,  and  sound-speed  variation.  Numerical  results 
indicate  that  some  current  structures  can  cause  large  variations  in  received 
intensity,  and  that  substantial  differences  can  occur  in  reciprocal  trans¬ 
missions.  Current  effects  on  intensity  may  be  quite  sensitive  to  the  sound- 
speed  distribution. 

INTRODUCTION 

Ocean  currents  cause  interesting  and  significant  effects  on  underwater 
sound.  For  example,  in  a  time-independent  ocean  environment,  currents  cause 
reciprocity  relations  to  fail.  Based  on  ray  theory,  work  has  been  done  by 
the  authors  and  others  to  model  the  influences  of  currents  on  acoustic 
transmissions.  It  has  been  shown,  for  instance,  that  certain  current  struc¬ 
tures  can  cause  large  fluctuations  in  total-field  intensity  and  per  ray 
phase.  These  results  are  limited  to  high-frequency  sound  transmissions.  To 
estimate  current-induced  effects  for  lower  frequencies,  a  full-wave  model 
should  be  used.  One  computational  model  is  the  parabolic  approximation, 
which  has  been  implemented  using  several  algorithms,  and  is  a  particularly 
attractive  method  for  efficiently  generating  transmission-loss  calculations. 

MODEL  FORMULATION 

In  a  complicated  medium  such  as  a  moving  ocean,  it  is  not  obvious  how 
known  parabolic  equations  should  be  modified  to  include  current  effects. 

For  this  reason,  we  systematically  reformulated  the  governing  time-dependent 
wave  equation,  starting  from  the  conservation  laws  and  state  relations 
governing  the  ocean  medium,  and  including  medium  motion  (Robertson  et  al., 
1985).  If  the  current  is  assumed  steady  but  non-uniform,  additional  terms 
occur  in  the  wave  equation.  These  terms  depend  on  the  current  gradient  and, 
as  will  be  discussed  below,  can  be  significant  in  subsequent  approximations. 
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Me  then  assume  that  the  sound  source  is  time -harmonic  and  that  the  propagat¬ 
ing  wave  is  outgoing,  thereby  obtaining  a  reduced  wave  equation.  After 
transforming  coordinate  systems,  we  invoke  the  far-field  approximation,  and 
proceed  to  generate  a  family  of  parabolic  equations.  Each  of  these  equa¬ 
tions  depends  on  the  relative  sires  of  three  dimensionless  parameters:  a 
Mach  number,  the  sound-speed  deviation,  and  a  shear  number.  The  last  pa¬ 
rameter  indicates  the  magnitude  of  current  gradient,  and  also  depends  on 
source  frequency.  For  example,  a  parabolic  approximation  appropriate  for 
an  isospeed  sound  channel,  through  which  flows  a  steady  depth-dependent 
current,  is 


2i + 
0Tr 


zz 


(1) 


where  Kq  is  a  wave  number,  u  is  current  speed,  and  is  an  envelope  of 
acoustic  pressure  in  the  far  field.  The  reference  sound  speed  cQ  in  this 
isospeed  case  is  equal  to  the  sound  speed  c.  Depth  is  indicated  by  z  and 
range  by  r.  It  can  be  shown  that  range-dependent  sound-speed  profiles  can 
be  incorporated  into  these  approximations  in  a  straightforward  way,  provided 
that  the  horizontal  gradient  of  sound  speed  is  not  large,  which  is  often  the 
case  in  many  important  ocean  regions.  For  example,  for  a  sound  channel  in 
which  the  sound-speed  deviation  is  of  the  same  order  of  magnitude  as  the 
Mach  number,  the  appropriate  equation  is 
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where  n(r,z)  =  cQ/c  is  the  index  of  refraction. 

In  Eqs.  (1)  and  (2),  the  term  which  depends  on  the  first  derivative  of 
current  may  or  may  not  be  retained,  depending  upon  the  size  of  the  shear- 
number  parameter.  In  the  event  that  it  is  kept,  the  parabolic  approxima¬ 
tions  are  not  in  a  "standard  form";  that  is,  they  cannot  be  solved  directly 
with  existing  numerical  implementations.  However,  it  is  possible  to  trans¬ 
form  this  family  into  related  parabolic  equations  which  are  in  a  standard 
form.  The  transformed  version  of  Bq.  (2)  is,  for  example: 

2i Kq <J>  +  $  +  k2  (n2-1)<f>  -  2k2  -  —  / d-2.u^ =  Or 

r  ZZ  0  0  \c  J  c0\dz2/  c2W/ 

where 


ij>(z,r)  =  e“Mu(z)  <f>(z,r )  •  (4) 

Note  the  appear ance  of  two  new  terms  in  Bq.  (3)  which  depend  on  the  square 
of  the  derivative  of  current  and  the  second  derivative  of  current.  The 
structure  of  Bq.  (3)  and  other  such  equations  suggests  the  use  of  an  effec¬ 
tive  sound-speed  profile  (ESSP),  which  includes  all  sound-speed  and  current- 
related  effects,  and  is  used  as  the  "actual"  sound  speed  for  numerical  solu¬ 
tion  of  the  equations.  For  example,  the  ESSP  corresponding  to  Bq.  (3)  is 


NUMERICAL  RESULTS 

To  solve  any  of  our  parabolic  equations  numerically,  we  elect  to  use 
the  IFD  model  developed  by  Lee  and  Botseas  (1982).  In  the  discussion  below, 
we  consider  several  sound-speed  profiles,  together  with  current  profiles,  as 
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Fig,  1.  Profiles  of  (a)  range-independent  sound  speed,  (b)  current, 
and  (c)  range-dependent  sound  speed, 

shown  in  Fig.  1.  The  zero-gradient  and  negative-gradient  sound-speed  pro¬ 
files  are  sometimes  range-independent  as  in  Fig.  1(a),  as  are  the  current 
profiles  in  Fig.  1(b).  The  surface  current  may  be  either  plus  or  minus  1 
m/s.  At  other  times,  we  employ  a  simple  range-dependent  sound-speed  pro¬ 
file,  Fig.  1(c),  for  which  the  horizontal  gradient  is  both  constant  and 
small.  Here,  isopleths  are  vertical  lines.  We  consider  first  the  result 
of  one  calculation  done  in  the  isospeed  channel,  with  surface  current  of 
magnitude  1  m/s.  The  source  frequency  is  200  Hz,  for  which  it  can  be  shown 
that  current-gradient  effects  are  negligible.  The  source  and  receiver 
depths  are  25  m.  In  Figs.  2-6,  the  bottom  acoustical  properties  are  the 
same  as  in  Robertson  et  al.  (1985).  Figure  2  shows  a  relative  intensity  in 
decibels  versus  range  for  three  cases:  no  current  present,  a  positive  cur¬ 
rent  in  the  source-receiver  direction,  and  a  negative  current  in  the  oppo¬ 
site  direction.  Several  important  current-related  effects  can  be  seen  in 
the  figure.  When  compared  to  the  solid  curve,  representing  intensity  in 
the  absence  of  any  current,  we  see  that  a  current  with  either  direction  can 
induce  substantial  variations  in  intensity.  For  example,  with  a  positive 
current  present,  variations  can  exceed  10  dB  over  certain  range  intervals, 
such  as  those  between  13  and  14  km.  Similar  behavior  is  seen  for  negative 
current.  Current  effects  in  both  cases  tend  to  increase  with  increasing 
range • 


Intensity  variations  are  highlighted  in  Fig.  3,  which  illustrates  dif¬ 
ference  in  relative  intensity  versus  range  for  three  cases:  no  current  and 
positive  current,  no  current  and  negative  current,  and  positive  and  negative 
currents.  Because  source  and  receiver  are  at  the  same  depth,  this  figure 
also  illustrates  one  type  of  effect  which  may  be  seen  in  reciprocal  trans¬ 
missions  (RTs).  The  intensity  difference  between  positive  and  negative  cur¬ 
rents,  indicated  in  Fig.  3  by  the  long-dashed  curve,  suggests  that  measure¬ 
ments  of  intensity  variation  between  reciprocal  source-receiver  pairs  may  be 
very  large.  Near  ranges  14  and  19  km,  this  difference  attains  a  magnitude 
of  nearly  20  dB.  At  other  range  intervals  the  difference  is  smaller,  but 
significant.  For  example,  between  10  and  12  km,  the  intensity  difference 
is  seen  to  generally  be  well  over  4  dB.  RT  differences  can  also  be  sig¬ 
nificant  in  range-dependent  channels.  Using  the  sound-speed  profile  in 
Fig.  1(c),  and  the  same  source  frequency,  source-receiver  depths,  and 
current  structure  as  above,  the  computed  intensity  difference  between  a 
source-receiver  pair  is  shown  in  Fig.  4.  Note  that  one  effect  of  the  range 


DIF r.  IN  NEl.  INTENSITY  (dB) 


Fig.  2.  Relative  intensity  versus  range  for  three  currents  in  an 
isospeed  channel • 


Fig.  3.  Difference  in  relative  intensity  versus  range,  from  Fig.  2. 

Oppositely-signed  current  effects  are  readily  compared  (long- 
dashed  curve ) • 
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Fig*  4.  Difference  in  relative  intensity  versus  range,  for  current 
reversal.  Reciprocal-transmission  difference  is  shown  for 
a  range-dependent  sound-speed  profile. 

variation  in  sound  speed  is  to  cause  the  relative  intensity  curve  to  shift 
toward  the  source,  when  compared  to  the  analogous  curve  in  Fig.  3.  Maximum 
intensity  values  have  been  altered,  also.  For  example,  at  a  range  14  km, 
the  peak  difference  is  over  30  dB. 

Another  way  to  visualize  the  effects  of  reversing  current  direction  is 
shown  in  Fig.  5.  Here,  differences  in  the  two  intensity  functions  are  plot¬ 
ted  as  level  curves  in  a  portion  of  the  depth-range  plane.  Ihe  difference 
is  intensity  for  a  positive  current  (in  the  source-receiver  direction)  minus 
intensity  for  a  negative  current,  as  in  Fig.  1(b).  The  source  frequency  and 
source  depth  are  again  200  Hz  and  25  m,  respectively.  Contour  intervals 
are  5  dB,  with  negative  differences  denoted  by  dotted  curves.  This  figure 
illustrates  the  intensity  differences  that  might  be  observed  in  a  channel 
with  tidal  effects.  At  ranges  larger  than  about  7  km,  bottom  attenuation 
has  stripped  away  most  higher  modes,  leaving  a  well-defined  pattern  of 
alternating  intensity  differences.  Regions  of  large  positive  difference 
occur  in  finger-like  patterns  which  alternately  emanate  from  the  channel 
surface  and  bottom.  Similar  structures  are  also  seen  for  negative  differ¬ 
ences.  In  this  example,  a  zone  of  very  small  differences  extends  in  range 
across  most  of  the  channel  at  roughly  mid— depth.  Regions  of  maximum  differ¬ 
ence  occur  regularly  above  and  below  this  zone.  In  contrast,  the  intensity- 
difference  pattern  in  a  negative  gradient  channel  is  noticeably  different. 
Figure  6  shows  the  analogous  level  curves  for  the  negative  gradient  profile 
of  Fig.  1(a).  Ihe  finger-like  structures  apparent  in  Fig.  5  have  in  some 
cases  blended  together  in  Fig.  6,  leaving  regions  of  high  intensity  differ¬ 
ence  located  at  many  mid-depth  points.  The  overall  pattern  is  more  compli¬ 
cated  than  the  one  present  in  the  isospeed  channel.  Consequently,  the 
intensity-difference  pattern  resulting  from  oppositely-signed  currents 
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Fig.  5.  Level  curves  of  intensity  difference  in  an  isospeed  channel, 
with  current  reversal. 


Fig.  6.  Level  curves  of  intensity  difference  for  a  negative  gradient 
channel  (see  Fig.  1(a)),  with  current  reversal. 
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appears  to  be  very  sensitive  to  these  types  of  changes  in  the  sound-speed 
profile. 

For  some  current  structures  and  lower  source  frequencies,  the  appropri¬ 
ate  transformed  parabolic  equation  will  include  new  terms  which  depend  on 
current  concavity  and  (possibly)  the  square  of  current  gradient.  One  type 
of  current  structure  which  may  require  additional  terms  is  shown  in  Fig. 
7(a).  At  the  surface,  the  current  speed  is  1  m/s,  and  it  decays  to  zero  at 
the  bottom.  Note  the  appearance  of  several  strong  shear  layers,  particular¬ 
ly  those  at  depths  35  and  60  m.  The  vertical  shear  structure  seen  here  can 
be  acoustically  significant  for  sufficiently  low  source  frequencies.  For 
example,  when  the  source  frequency  is  30  Hz,  the  ESSP  is  similiar  to  that 
given  by  Bg.  (5)  and  is  depicted  in  Fig.  7(b).  In  this  example,  concavity 
effects  are  significant,  but  shear  effects  can  be  neglected.  Note  that 
current  concavity  dominates  the  behavior  of  the  ESSP.  lhe  current  shear 
structure  has  introduced  large  rapid  variations,  one  of  which,  at  the  depth 
60  m,  approaches  20  m/s  in  Fig.  7(b).  For  higher  source  frequencies,  the 
magnitude  of  the  variations  decreases,  yet  may  still  be  significant.  We 
anticipate  that  this  current  structure  can  cause  interesting  acoustical 
effects . 

In  Fig.  8,  we  see  one  result  of  computations  done  with  the  current 
structure  shown  in  Fig.  7(a).  In  order  to  observe  concavity  (or  second 
derivative)  effects,  we  solved  the  relevant  parabolic  equation,  first  with 
concavity  included  and  then  with  concavity  omitted.  The  root-mean-square 
difference  of  the  intensities  in  the  two  cases,  called  J,  was  then  calcu¬ 
lated  with  a  range  averaging.  Source  and  receiver  depths  are  25  m,  and  the 
source  frequency  is  100  Hz.  The  results  for  three  different  bottom  types 
are  shown.  The  rigid  bottom  is  perfectly  reflecting,  while  beneath  the 
water  column  for  the  hard  and  soft  bottoms  is  a  second  fluid  layer  with  dif¬ 
ferent  sound  speed  and  density,  lhe  hard  bottom  has  larger  discontinuities 
in  these  quantities  than  the  soft  bottom.  For  both  hard  and  soft  bottoms, 
a  small  amount  of  volume  attenuation  was  introduced.  As  the  bottom  changes 
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Fig.  7.  Profiles  of  (a)  a  current  with  high  shear,  and  (b)  the 
effective  sound  speed. 
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Fig.  8.  RMS  difference  in  intensity  versus  range  for  the  current  of 
Fig.  7(a)  and  for  three  bottom  types. 

from  rigid  to  hard  to  soft,  note  that  the  overall  values  of  J  tend  to  de¬ 
crease.  However,  the  peak  values  may  actually  increase  substantially.  For 
example,  at  45  km  the  soft  bottom  has  a  peak  which  is  about  4  dB  larger  than 
its  overall  value.  Furthermore,  the  curves  are  smoother  for  both  the  hard 
and  soft  bottoms,  since  they  attenuate  higher  modes  more  rapidly  than  the 
rigid  bottom.  These  observations  illustrate  the  strong  dependence  of  con¬ 
cavity  effects  on  bottom  influences  in  underwater  sound  transmissions. 

SUMMARY 

We  discuss  a  family  of  parabolic  approximations,  valid  for  depth-  and 
range-dependent  sound-speed  profiles,  which  include  effects  caused  by  steady 
depth-dependent  currents.  These  approximations  permit  examination  of  inten¬ 
sity  effects  caused  by  currents  for  frequencies  and  environments  where  other 
models  may  not  be  valid  or  convenient.  Using  a  standard  numerical  implemen¬ 
tation,  we  present  the  results  of  computations  for  several  current  and  sound- 
speed  structures.  They  suggest  that  currents  can  cause  significant  intensi¬ 
ty  variations,  principally  by  altering  the  effective  sound-speed  profile. 
Intensity  differences  arising  from  reciprocal  transmissions  are  shown  to  be 
especially  large.  Also,  current  effects  on  intensity  can  be  very  sensitive 
to  small  changes  in  sound  speed.  Finally,  the  presence  of  current  fine 
structure  can  introduce  additional  fluctuations  in  intensity  predictions. 
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Defense  Advsnced  Research  Projects  Agency 
1400  Wilson  Boulevard 
Arlington,  Virginia  22209 

Attni  Capt.  V.  Si  None  1 

ARPA  Research  Center 
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CoNsnding  Officer 
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U.S.  Naval  Academy 
Annapolis,  Maryland  21402 
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Applied  Research  Laboratories 
University  of  Texas  at  Austin 
P.0.  Box  0029 
10000  FM  Rood  1325 
Austin,  Toxoo  78712 
Attm  Dr.  Loyd  Hampton 
Dr.  Chariot  Wood 

Atlantic  Oceanographic  and 
Meteorological  Laboratoriaa 
15  Rickenbecker  Causaway 
Hi  ami «  Florida  33149 
Attn*  Dr.  John  Proni 

Ball  Telaphone  Laboratoriaa 
1  Whippeny  Road 
Whippeny,  New  Jersey  07981 
Attn:  Dr.  Bruce  Bogart 
Dr.  Pater  Hiroch 

Bolt,  Beranek,  and  Nemsan,  Inc. 

50  Moulton  Street 

Cambridge,  Massachusetts  02238 

Attn*  Dr.  K.  L.  Chendiramani 

Chase,  Inc. 

14  Pinckney  Street 
Boston,  Massachusetts  02114 
Attn:  Dr.  David  Chase 

Dr.  David  Middleton 
127  Cast  91st  Street 
New  York,  New  York  10028 

Duke  University 

Department  of  Electrical  Engineering 
Durham,  North  Carolina  27706 
Attn:  Dr.  Loren  Nolte 

General  Electric  Company 
Heavy  Military  Electronic  Systems 
Syracuse,  New  York  13201 
Attn:  *.  Don  Winfield 
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Schenectady,  New  York  12301 
Attn*  Dr.  Thomas  G.  Kincaid 

Gould,  Incorporated 
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6711  Baymeadow  Drive 
Glen  Bumie,  Maryland  21061 
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Hughes  Aircraft  Company 
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321  Northland  Ava. 

P.0.  Box  3618 

Rochester  i  New  York  14610 

Institute  for  Acoustical  Research 
Miami  Division  for  the  Palisades 
Geophysical  Institute 
615  South  West  2nd  Avenue 
Miami,  Florida  33130 
Attn*  i  Hr.  M.  Kronengold 
Dr.  J.  Clark 

Institute  of  Geophysics  and  Planetary 
Physics 

Scripps  Institute  of  Oceanography 
University  of  California 
La  Jolla,  California  92093 
Attn:  Dr.  W.  Munk 

Mr.  J.  Spiesberger 
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205  South  Whiting  Street 
Suita  607 

Alexandria,  Virginia  22304 
Attn:  Dr.  S*  Ada 
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Massachusetts  Institute  of  Technology 
Acoustics  and  Vibration  Laboratory 
70  Maasachusatta  Avenue 
Room  5-222 

Cambridge,  Massachusetts  02139 
Attn:  Professor  Patrick  Leehey 

Palisades  Sofar  Station 
Beneuda  Division  of  Palisades 
Geophysical  Institute 
FP0  New  York  09560 
Attn:  tfr.  Carl  Hartdegen 

Polar  Research  Laboratory 
123  Santa  Barbara  Avenue 
Santa  Barbara,  California  93101 
Attn:  Mr.  Beaumont  Buck 

Research  Triangle  Institute 
Research  Triangle  Park 
Durham,  North  Carolina  27709 
At  try  Dr.  S.  Huffman 

Rensselaer  Polytechnic  Institute 

Troy,  New  York  12181 

Attn:  Dr.  Melvin  J.  Jacobson 

Science  Applications,  Inc. 

8400  Westpark  Drive 
McLean,  Virginia  22102 
Attn:  Dr.  P.  Tatro 

S.D.P.  Inc. 

15250  Ventura  Boulevard 
Suita  518 

Sherman  Oaks,  California  91403 
Attn:  Dr.  M.  A.  Basin 
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Tun  Instnaients,  Inc. 

13500  North  Central  Expressway 
Dallas,  Texas  75231 
Attni  tt.  Charles  Black 

Underwater  Systeas,  Inc. 

6121  Georgia  Avenue 

Silver  Spring,  Maryland  20910 

AtLns  Dr.  M.  Weinstein 

University  of  Mi  sal 
Rosenatlel  School  of  Marine  and 
Ateospheric  Sciences 
4600  Rickenbacker  Causeway 
hiaai,  Florida  33149 
Attn:  Dr.  H.  DeFerrari 

University  of  Michigan 
Departaent  of  Aerospace  engineering, 
North  Caw>us 

Ann  Arbor,  Michigan  4B109 
Attn:  Dr.  W.  W.  Wilaarth 

University  of  Michigan 
Cooley  Electronics  Laboratory 
Ann  Arbor,  Michigan  48105 
Attn:  Dr.  T.  C.  Birdsall 

University  of  Rhode  Island 
Departaent  of  Electrical  Engineering 
Wakefield,  Rhode  Island  02881 
Attn:  Dr.  Donald  Tufts 

Woods  Hole  Oceanographic  Institution 
Woods  Hole,  Massachusetts  02543 
Attn:  Dr.  Paul  McElroy 
Dr.  R.  Spindel 
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